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I.  INTRODUCTION 

The  Solid  Fuel  Ramjet  (SFRJ)  family  of  projectiles  is  a  new  technology  in  the  Field  of  ballis¬ 
tics  that  features  a  capability  to  modulate  the  acceleration  in  flight.  At  present,  SFRJ  training 
rounds  are  being  developed  that  simulate  low  drug  at  the  initial  phase  of  the  flight  and  produce 
high  deceleration  past  the  burning  time.  The  SFRJ  famity  of  projectiles  consists  of  spin- 
stabilized  and  fin-stabilized  configurations.  Both  types  have  Pitot  type  inlets.  The  fln-stabilized 
configurations  have  circumferential  ridges  or  grooves  that  transfer  the  launch  accelerating  force. 

The  Missile  DATCOM1  and  the  NSWC  Aeroprediction2  codes,  that  are  available  at  BRL.  are 
capable  of  analyzing  the  aerodynamic  characteristics  of  body  alone  and  ftn-body  combinations. 
However,  neither  code  includes  an  option  of  an  inlet  in  the  body  or  takes  into  account  the  aero¬ 
dynamic  effects  associated  with  the  roughness  elements.  Since  the  codes  possess  quick  response 
capability,  a  hybrid  scheme  of  analysis  is  devised.  The  available  codes  axe  used  for  basic  analy¬ 
sis  of  model  configurations.  Then,  corrections  for  the  two  effects  that  were  mentioned  above,  are 
applied. 

Selection  cf  an  appropriate  prediction  code  is  considered  first.  Then,  the  additional  analysis  is 
described.  Finally,  the  application  of  the  proposed  scheme  to  the  various  SFRJ  configurations  is 
presented.  The  predictions  are  compared  with  available  test  data  from  wind-tunnels  and  ballistic 
ranges. 


II.  SELECTION  OF  A  DESIGN  CODE 

The  Missile  DATCOM  and  the  NSWC-Aeroptediction  codes  are  capable  of  estimating  the 
external  axial  force,  the  normal-force  and  the  static  stability  of  the  pertinent  configurations.  For 
the  analysis  of  body  alone,  in  the  supersonic  domain,  both  codes  offer  the  same  analytical 
options:  Van  Dyke’s3  hybrid  theory  and  the  second  order  shock  expansion4  method  (SOSE). 
Indeed,  the  two  codes  give  practically  the  same  answers  when  applied  to  the  same  body  con¬ 
figuration.  Only  the  SOSE  can  cover  the  whole  supersonic  Mach  number  range.  Comparison  of 
prediction  by  this  method  to  empirical  data  bases  showed  very  good  agreement  in  normal-force 
curve  slope  and  center-of-pressure.  They  will  be  discussed  in  Chapter  IV,  with  regard  10  actual 
configurations. 

Both  codes  use  the  supersonic  lifting  surface  method,  when  applicable,  for  the  fin  alone 
aerodynamic  characteristics. 

The  main  difference  between  the  codes,  for  the  present  application,  is  in  the  body-fin  and 
fin-body  influence  coefficient.  For  fins  mounted  entirely  on  a  body,  both  codes  use  Nielsen's5 
influence  coefficients.  For  fins  that  over-hang  the  base  of  the  body,  as  in  the  case  of  the  60mm 
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SFRJ-F1NNER  projectile  that  will  be  described  in  Chapter  IV,  Section  3,  the  DATCOM  code 
uses  the  same  influence  coefficients.  Only  the  NSWC  code  applies  a  correction,  bused  on 
empirical  data,  to  account  for  the  reduction  in  interference  normal-force  associated  with  the 
over-hang.  As  will  be  elaborated  below,  the  correction  is  not  small  so  that  the  NSWC  code  has 
an  advantage.  Also,  base  drag  coefficients  predicted  by  this  code  are  in  good  agreement  with  the 
flight  data  of  Stoney* 


111.  POSTPROCESSING 


Since  the  NSWC  code  does  not  accept  inlets,  the  tangent  ogive  cowl,  in  ^e  computational 
model,  is  extended  to  form  a  closed,  pointed  nose.  Assuming  that  slender  body  theory  is  valid 
for  the  nose  extension,  no  correction  is  made  for  the  normal-force  slope  as  the  contribution, 


(1) 


is  equal  to  that  predicted  by  momentum  change  of  the  ingested  stream-tube.  Furthermore,  it  is 
assumed  that  the  contribution  of  the  actual  cowl  to  is  the  same  as  with  nose  extension.  The 

center-of-pressure,  however,  is  shifted  backward.  The  amount  of  shift  depends  on  the  structure 
of  the  shock  and  expansion  waves  within  the  inlet.  Since  this  is  not  known  at  present,  a  conser¬ 
vative  assumption  is  used,  namely  that  the  above  contribution  acts  at  the  inlet  plane.  Approxi¬ 
mating  the  foremost  part  of  the  extended  nose  by  a  cone,  having  a  length  l„  the  change  in  pitch¬ 
ing  moment  coefficient  is: 

ACMou  *  ZJ  CVcu  ’ 


The  associated  shift  in  center-of-pressure  is: 


(2) 


(3) 


The  roughness  elements,  required  by  sabot  launch,  generate  additional  drag  that  depends  on 
geometry,  Mach  and  Reynolds  numbers  and  flight  conditions.  The  additional  drag  reduces  the 
dynamic  pressure  of  the  flow  approaching  the  fins.  It  is  expected  that  the  reduced  dynamic 
pressure  will  decrease  the  efficiency  of  the  fins  as  stabilizers. 
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The  full  treatment  of  the  problem  is  a  topic  for  computational  fluid  dynamics.  Here,  an 
engineering  approximation  is  proposed,  for  use  in  conjunction  with  the  design  code. 

The  drag  coefficient  of  the  roughness  elements  can  be  estimated  based  on  correlations  of 
experimental  data  such  as  by  Hoemer7  and  Young  and  Patterson;8 

*  4S  CD'  W 


where  the  superscript  s  refers  to  a  series  of  protuberances  or  grooves,  considering  the  interfer¬ 
ence  between  the  elements. 

The  additional  axial-force  appears  as  an  increase  in  boundary-layer  momentum  thickness 

*  J  *  J  *i>J  •  l» 


The  thickening  of  the  boundary  layer  also  causes  an  increase  in  the  displacement  thickness, 
A8*.  In  the  case  of  the  60mm  SFRJ-FINNER,  that  will  be  discussed  in  the  next  chapter,  the 
thickness  of  the.  boundary-layer  ahead  of  the  roughness  is  very  small  compared  to  the  increase 
due  to  the  threads.  It  is  assumed  that  the  shape  factor  over  the  rough  surface  is  the  same  as  that 
over  a  smooth  one.  The  shape  factor,  H  «  A5*/A8**  .depends  on  the  velocity  and  density 
distribution.  A  common  assumption  is  the  validity  of  the  1/7  power  law  at  supersonic  Mach 
numbers  and  adiabatic  flow.  This  assumption  will  be  employed  in  the  analysis  that  follows. 

The  displacement  thickness,  with  roughness  elements,  is  not  small  compared  to  the  diameter 
of  the  body  or  the  exposed  span  of  the  fins.  Two  aerodynamic  effects  are  produced  by  the 
thickened  boundary  layer.  The  first  one  is  the  flaring  of  the  center-body  to  an  effective  diameter 
D»  »  (D  +  2A8*).  The  contribution  of  the  effective  flare  to  normal-force  curve  slope,  as  pre¬ 
dicted  by  slender  body  theory  is; 

=  {(&)  - ']  • 

In  the  present  case,  only  the  actual  cylindrical  body  carries  normal  loads.  Hence,  the  above 
expression  will  be  multiplied  by  the  ratio  of  the  cylinder  projected  area  to  that  of  tne  effective 
flare.  Introducing  the  effective  flare  diameter,  the  additional  contribution  of  the  body  due  to 
roughness  is: 
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^or  =  21 


r*  /■  +  \2 

D  +  2A5  A 


-  I 


D  _  0A5 


n 


The  second  effect  is  the  reduction  of  the  contribution  of  the  fins  to  the  nonml-force  due  to 
masking  by  the  body  boundary-layer.  It  is  assumed  that  the  inner  sections  of  the  fins,  with  width 
equal  to  A5  , income  inefficient.  The  mask  area  on  a  pair  of  fins  is  Sm  ~  2 Cr  Ao  .  Thus,  the 
contribution  of  a  fin  set  to  the  normal-force  and  the  pitching-moment  is  reduced  by  a  factor 


S-S 


n  = 


m 


(6) 


s  1-2 


CR  A5 


The  corrected  urmal-force  curve  slope  and  the  center-of-pressure  location  are 


CNa*  *  cNa*  +  nC^a/  +  ACNar 


CMa  =  cMah  *  n cMaf  +  *cMai  +  acMoj 


+  r  + 
Xcp  _  CWa 

D  "ch  ■ 
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IV.  APPLICATION 


l.  40mm  SFRJ  PROJECTILE 

The  external  dimensions  of  the  spin  stabilized  40mm  SFRJ  configuration  are  shown  in  Figure 
1.  The  truncated  ogival  nose  has  a  fineness  ratio  of  1.8  and  the  cylindrical  center  body  a  fineness 
of  ratio  2.4.  The  inlet  diameter  ratio  is  D  JD  =  0.5 1 . 

The  computed  normal-force  curve  ,sh  pe  and  center-of-pressure  art  shown  in  Figure  2a  and 
Figure  2b,  respectively.  The  Van  Dyke  method  predicts  a  smaller  lift  carry  over  than  that 
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predicted  by  the  SOSE  method.  It  shows  as  a  smaller  normal-force  slope  and  a  more  forward 
center-of- pressure.  The  correction  for  the  inlet  opening,  relative  to  the  computational  model  that 
has  an  extended  nose,  amounts  to  a  buckward  shift  of  center-of- pressure  of  about  0.045  calibers. 

Comparison  of  the  prediction  with  range  flight  data,9  and  wind-tunnel  data,10  is  also  presented 
in  Figure  2.  The  measured  slope  of  the  normal-force  curve  is  18-23%  higher  than  that  calculated. 
Thus,  the  experimental  Aircraft  DATCOM11  data  base  is  used  for  further  validation  of  the 
analysis.  This  empirically  based  source,  substantiates  the  results  of  the  SOSE  method,  for 
extended  nose  body,  down  to  a  Mach  number  of  1 .5.  It  is  concluded  that  the  cowl  produces  more 
normal-force  with  an  open  inlet  than  with  an  extended  nose.  It  should  be  noted  that  the  range 
data  is  an  average  of  15  tests  with  large  scatter  and  an  estimated  error  in  excess  of  30%. 

The  data  obtained  in  free  flight  shows  an  average  center-of-pressure  location  which  is  0.08 
calibers  more  aft  than  that  predicted.  The  wind-tunnel  data  are  in  near  perfect  agreement  with 
the  analysis  at  low  Mach  numbers  and  gives  up  to  0.24  calibers  in  a  more  aft  location  at  the  high 
end  of  the  Mach  number  range. 

2.  75mm  SFRJ  PROJECTILE 


The  75mm  SFRJ  is  a  spin  stabilized  projectile  having  a  nose  Fineness  ratio  of  2.0  and  a  center 
body  fineness  ratio  of  1.6.  The  inlet  diameter  ratio  is  0.474  (also  see  Figure  3). 
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The  analytically  computed  normal-force  curve  slope  and  center-of-pressure  are  presented  in 
Figure  4a  and  Figure  4b,  respectively.  Here,  the  correction  to  the  pitching-moment  coefficient 
shifts  the  center-of-pressure  backward  about  0.04  calibers.  As  in  the  case  of  the  40mm  projec¬ 
tile,  no  boundary  layer  displacement  thickness  corrections  were  applied. 


The  agreement  between  the  predicted  static  aerodynamic  coefficients  and  those  obtained  at  a 
Mach  number  4.44  in  ballistic  range  tests12  is  very  good,  as  is  apparent  in  Figure  4.  The  ex- 
j  perimentally  obtained  normal-force  curve  slope  is  only  3.7%  higher  than  that  predicted  and  the 

center-of-pressure  location  is  0.1  calibers  more  forward. 

|  3.  60mm  SFRJ-FINNER  PROJECTILE 

L  The  60mm  SFRJ-FINNER  projectile,  shown  in  Figure  5,  is  sabot  launched.  The  truncated 

ogival  nose  has  a  fineness  ratio  1.63  and  the  inlet  diameter  ratio  is  DJD  =  0.5745.  The  center 
body  has  fineness  ratio  5.60  and  features  four  pairs  of  buttress  threads.  Two  sets  of  fins  were 

(used,  as  shown  in  Figure  5a  and  Figure  5b.  The  original  one  has  aspect  ratio  of  0.92  and  the 
extended  one  aspect  ratio  of  0.528.  Both  sets  have  the  same  exposed  span  which  is  0.9745  body 
diameters. 
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As  mentioned  in  Chapter  II,  the  main  difference  between  the  two  candidate  codes  is  in  the 
influence  coefficient.  This  difference  is  illustrated  in  Figure  6  for  the  present  configuration  with 
original  fins  at  Mach  numbers  of  2  and  4.  The  fin  influence  coefficient,  is  assumed 

independent  of  the  fin  over-hang  by  the  DATCOM  code.  On  the  other  hand,  the  NSWC  code 
predicts  a  decrease  in  the  coefficient  as  over-hang  ratio  increases.  The  body  influence  coeffi¬ 
cient,  predicted  by  the  NSWC  code  is  somewhat  smaller  than  that  predicted  by  the 

DATCOM  code.  Both  codes  assume  that  this  influence  coefficient  diminishes  as  the  square  of 
one  minus  the  over-hang  ratio.  The  sum  of  the  two  influence  coefficients  predicted  by  the 
NSWC  code  is  about  15%  lower  than  that  predicted  by  the  DATCOM  code. 


A  step  in  the  validation  of  the  code  was  a  comparison  of  computed  body  alone  characteristics 
with  data  bases.  The  Aircraft  DATCOM,11  used  in  Section  1,  and  the  ESDU13  data  bases  were 
used.  The  comparison,  for  the  extended  nose  body,  is  shown  in  Figure  7.  The  analysis  is  in 
excellent  agreement  with  the  data  bases  at  Mach  numbers  between  2.0  and  3.0.  The  code  under¬ 
predicts  the  normal-force  curve  slope  by  6.7%  at  Mach  number  of  1.5  and  by  8.3%  at  5.0.  The 
predicted  center-of- pressure  location  is  0.1  calibers  more  forward  than  that  predicted  by  the  data 
base  at  a  Mach  number  of  1.5  and  0.14  calibers  at  5.0.  To  summarize  this  comparison,  as  in  the 
case  of  the  40mm  SFRJ,  the  data  bases  corroborate  the  prediction  by  the  SOSE  method  down  to  a 
Mach  number  of  1 .5. 

The  drag  coefficient  of  the  buttress  threads  was  estimated  using  the  empirical  data  compiled 
by  Young  and  Patterson.8  It  is  assumed  that  the  drag  produced  by  each  thread  is  a  sum  of  a 
chamfered  forward  facing  step  and  that  of  a  backward  facing  one.  Each  of  these  contributions 
depend  upon  flight  conditions  and  the  smooth  wall  skin  friction.  Since  data  is  limited,  the 
dependence  on  Mach  number  is  based  on  that  of  a  straight  ridge.  For  this  protuberance,  it  was 
found,  from  Figure  4.37  of  Reference  8,  for  1.4  £  M  <,  2.8  that  the  drag  coefficient  ratio  for  a 
constant  Reynolds  number  is: 


WM)  (m  X* 

2J 


(7) 


Since  no  data  were  found  for  M  >  2.8,  'his  expression  was  used  up  to  M  =  5.0.  The  drag  coeffi¬ 
cient  curve  of  a  single  thread,  based  on  frontal  area,  Cp»  ,  was  evaluated  for  two  flow  conditions. 

e 

One  is  for  full-scale  flight  at  sea  level,  assuming  that  the  surface  reaches  recovery  wall  tempera¬ 
ture;  the  other  is  for  NASA  Langley  Research  Center  (LaRC),  unitary  wind-tunnel,  where  tests, 
which  will  be  described  in  the  following  chapter,  were  performed.  The  estimates  are  presented  in 
Figure  8  and  show  that  the  drag  coefficient  of  a  single  thread  under  the  LaRC  wind-tunnel  flow 
conditions  is  much  smaller  than  that  iii  free  flight.  This  difference  is  caused  by  the  large  differ¬ 
ence  in  Reynolds  number  between  the  two  cases. 
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The  contributions  of  the  eight  threads  in  series  is  evaluated  using  the  empirical  data  of  Lacey, 
reported  in  Figure  4.46  of  Reference  8.  It  is  estimated  that  for  the  present  average  ridge  density 
and  surface  length  to  ridge  height  ratio,  the  drag  of  the  series  is  3.74  times  that  of  a  single 
isolated  thread.  It  is  assumed  that  this  factor  is  independent  of  Mach  number.  The  additional 
axial-force  coefficient,  due  to  the  series  is  evaluated  using  Equation  (4)  and  presented  in  Figure  9 
for  the  two  flow  conditions  described  above. 

As  mentioned  above,  the  boundary  layer  shape  factor  is  based  on  the  assumption  of  the 
validity  of  the  power  law  at  supersonic  Mach  numbers  and  adiabatic  conditions,  and  that  it  is 
unaffected  by  surface  roughness.  These  assumptions  were  validated  by  comparing  experimen¬ 
tally  obtained  data  by  Fenter14  (sand  roughness,  1.62  <  M  <  2.70)  and  Voisinet15  (screens,  M  = 
2.95)  to  the  analytical  results  of  Shapiro.16  The  agreement  is  good  and  the  differences  between 
smooth  and  rough  wall  cases  are  small,  as  shown  in  Figure  10. 

The  momentum  thickness,  associated  with  the  additive  drag  due  to  the  buttress  threads, 
decreases  as  Mach  number  increases,  as  shown  in  Figure  11.  For  the  free-flight  case,  the  dis¬ 
placement  thickness  increases  with  Mach  number  and  reaches  a  value  of  0.05  bodies  diameters  at 

M  *  5.0.  For  the  LaRC  flow  conditions,  8*  decreases  as  the  Mach  number  increases  and  its 
value  at  M  *  5.0  is  0.017  body  diameters.  The  associated  fin  efficiencies  is  shown  in  Figure  12. 
For  M  5.0,  which  is  the  design  Mach  number,  the  fin  efficiencies  at  flight  conditions  are  0.83 
and  0.87  for  the  original  and  extended  fin  configurations,  respectively.  The  difference  between 
the  two  cases  is  due  to  the  variance  in  their  aspect-ratios  and  the  taper  ratios.  In  the  LaRC  flow 
conditions,  the  fin  efficiency  is  very  close  to  1.0,  reflecting  the  small  axial-force  due  to  the 
buttress  threads. 

Figure  13  shows  the  predicted  aerodynamic  characteristics  of  the  configuration  having  the 

original  fins  in  free-flight  conditions.  The  code  predicts  instability  at  Mach  numbers  larger  than 

2.0.  A  range  test  at  a  Mach  number  of  4.25  showed  that  this  configuration  is  indeed  unstable. 

The  broken  lines  between  Mach  numbers  1.5  and  2.0  indicate  computational  transition  from  the 

use  of  a  transonic  method  to  the  use  of  supersonic  lifting  surface  method  in  the  evaluation  of  fin 

alone  characteristics.  In  the  present  case,  the  code  uses  the  value  C ^  =  irAR/2,  as  predicted 

ot 

by  slender  wir.'i  theory,  for  Mach  numbers  for  which  AR(3  <  2.0.  The  jump  between  this  value 
and  that  predicted  by  supersonic  lifting  surface  method  is  not  physical.  It  is  expected  that  the 
actual  curves  will  be  more  smooth  that.  those  predicted  by  the  code. 

Results  of  the  analysis  for  the  configuration  having  the  extended  fin,  under  free-flight  condi¬ 
tions,  are  shown  in  Figure  14.  The  normal-force  curve  slope  (Figure  14a)  with  buttress  threads  is 
almost  equal  to  that  of  the  smooth  body  for  Mach  numbers  smaller  than  3.5  and  slightly  higher  at 
Mach  numbers  larger  than  4.0.  The  fprward  shift  of  center-of-pressure  (Figure  14b)  due  to  the 
threads  reaches  a  value  of  0.20  calibers  at  a  Mach  number  of  5.0.  The  analysis  predicts  a  very 
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small  margin  of  stability  for  a  smooth  body  and  instability  at  Mach  numbers  larger  than  4.3  for 
the  rough  one. 

Results  of  one  range  flight  test  were  available  at  the  time  of  preparation  of  this  report.  The 
measured  center-of-pressure  falls  between  the  two  predicted  curves.  The  measured  normal-force 
curve  is  nine  percent  higher  than  that  computed. 


V.  WIND-TUNNEL  TEST  OF  THE  60mm  SFRJ-FINNER  PROJECTILE 


1.  OBJECTIVES 


|  The  objectives  of  the  two  sets  of  wind-tunnel  tests,  that  will  be  described  in  this  chapter  are: 

I 

I 

1  1.  To  determine  the  margin  of  static  stability  of  the  60mm  SFRJ-FINNER  projectile  as  a 

!  function  of  Mach  number. 


2.  To  determine  the  effects  of  the  buttress  threads  on  SFRJ  projectile  stability. 

3.  To  determine  the  effect  of  ingested  inlet  flow  on  normal-force  and  stability. 


4.  To  calibrate  current  design  codes  in  terms  of  their  ability  to  predict  SFRJ  projectile 
aerodynamic  characteristics. 

Other  factors  considered,  such  as  grooves  versus  buttress  threads,  fin  over-hang  and  fin  size  were 
tested  but  will  not  be  considered  heie. 

2.  EXPERIMENTS 

The  NASA  Langley  Research  Center  (LaRC)  Unitary  wind-tunnel  is  an  asymmetric,  4  ft.  by  4 
ft.  test  section,  facility.  The  test  plan  by  Yalamanchili17  called  for  force  and  moment  measure¬ 
ments  at  five  Mach  numbers  between  3.0  and  4.6  with  stagnation  conditions  adjusted  to  give  a 
Reynolds  number  per  foot  of  2  x  106.  The  model  was  basically  a  full  scale  60  mm  SFRJ  with 
buttress  threads.  A  second  model  was  tested  with  the  space  between  the  threads  filled  so  as  to 
produce  i  smooth  continuous  body  shape.  Tests  were  also  performed  with  an  extended  nose 
attached,  eliminating  the  flow  through  the  model.  A  six  component  internal  strain  gauge  balance 
was  used  with  sting  mounting  of  the  model.  When  operated  with  an  open  inlet,  the  flow  was 
believed  swallowed  although  Schlieren  photography  failed  and  was  not  available. 

The  US  Army  Chemical  Research,  Development  and  Engineering  Center  (CRDEC)  tests  by 
Olsen  et.  al."  were  performed  in  the  Aerodynamics  Research  and  Concepts  Assistance  Pranch 
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supersonic  wind-tunnel.  This  tunnel  is  also  an  asymmetric,  blow  down  tunnel  with  a  6  in.  by  6 
in.  test  section.  The  tunnel  was  operated  at  lower  but  overlapping  Mach  numbers  (2.00  to  3.75) 
with  respect  to  the  LaRC  experiment.  The  air  supply  conditions  involved  higher  stagnation 
pressures  so  that  the  Reynolds  number  per  foot  were  five  times  the  Langley  value.  Because  the 
CRDEC  tunnel  is  considerably  smaller,  a  0.3175  scale  model  of  the  60mm  SFRJ  configuration 
was  employed.  The  model  was  also  sting  mounted  with  an  internal  strain  gauge  balance  for 
measurement  of  normal-force  and  pitching  moment. 

Both  facilities  reduced  their  measurements  to  normal-force  and  pitching-moment  coefficients 
versus  angle  of  attack.  It  is  assumed  that  sting  deflection  and  instrument  errors,  if  any,  have  been 
accounted  for  in  the  reduced  data. 

3.  DATA  ANALYSIS 

The  vs  a  and  vs  a  data  were  analyzed  by  fitting  an  expression  of  the  form 

y  =  a  +  bet  +  calrd  . 

This  equation  was  selected  because  it  is  the  simplest,  low  order,  formula  which  displays  the 
general  form  observed  and  expected  in  the  data.  That  is,  ideally  (perfectly  symmetric  model  in 
uniform  flow)  the  forces  and  moments  are  expected  to  be  odd  functions  of  a.  However,  some 
non-ideal  behavior  is  observed  in  that  the  forces  and  moments  do  not,  as  a  rule,  go  to  zero  at  zero 
angle  of  attack.  The  values  of  a,  b  and  c  are  determined  by  a  least  square  technique.  The  best  fit 
value  of  b  is  interpreted  as  the  slope  of  the  normal-force  or  pitching-moment  coefficient  curve, 
dC^/da  or  dC^j/8a ,  depending  on  whether  y  =  Cj^  or  C^|.  The  center  of  pressure  is  then 

calculated  from  these  slopes: 

Xcp  Cwa 

D  ~  CV 

''a 


where  the  X  is  measured  from  the  reference  point  in  the  definition  of  CM.  The  LaRC  data 

referenced  the  moments  to  a  point  0.5358L  (L  =  19.03  inches  or  the  length  of  the  closed  inlet 
configuration)  from  the  closed  inlet  nose.  The  CRDEC  reference  for  moments  is  the  model  base. 
In  the  results  reported  here  Xcp  is  referenced  to  the  lip  of  the  open  cowl  configuration  even  for  the 
closed  inlet  case. 


4.  ACCURACY 

The  accuracy  of  the  Langley  balance  is  cited  as  ±0.25  percent  of  full  scale  which,  for  each 
component,  is  ±0.25  Ibf  or  ±0.5  lbfiin.  for  normal-force  and  pitching  moment,  respectively. 
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These  values  can  be  rt presented  as  percem  uncertainties  in  N  and  M.  For  example,  at  »  3.0 
and  a  =  6°,  N  =  8.2  lbf  ±3%  and  M  *  14.4  lbf-in.  ±3.5%. 

An  alternate  evaluation  of  the  uncertainties  can  be  derived  from  the  statistics  of  the  curve  Fits. 
The  analysis  can  be  extended  to  all  the  aerodynamic  characteristics  (C^,  ^V/a  ant*  *cp)-  From 

each  curve  fit  an  estimated  standard  deviation  is  calculated  from: 

n 

o  =  {X  [yj  -  (a  +  betj  +  ccxj|ctj|)|2/(n  -  3)}1/2 
3  i=I 

A  plot  of  a  sample  versus  a  and  versus  a  are  given  in  Figure  15.  These  results  are 

for  M oo  *  3.00,  and  a  closed  inlet  configuration  with  buttress  threads.  The  curve  fits  are  shown 

as  solid  or  dashed  lines  for  the  LaRC  or  CRDEC  data,  respectively.  The  estimated  standard 
deviations  for  this  typical  case  are: 


or 

CM 

LaRC 

0.0041 

0.0045 

CRDEC 

0.013 

0.039 

Based  on  the  averages  of  the  estimated  uncertainties  of  the  LaRC  and  CRDEC  results.  The 
CRDEC  data  have  three  times  more  scatter  and  corresponding  uncertainty.  This  difference  can 
be  directly  attributed  to  the  smaller  scale  of  the  model  in  the  CRDEC  experiment. 

5.  STABILITY  OF  THE  60mm  SFRJ-F1NNER  PROJECTILE 

Figure  16  shows  the  experimentally  obtained  center-of-pressure  location  of  the  60mm  SFRJ. 
Only  the  actual  configuration,  namely  the  one  with  open  inlet  and  buttress  threads,  is  included. 
In  general,  the  LaRC  and  CRDEC  results  are  in  good  agreement  despite  the  significantly  larger 
uncertainty  in  the  CRDEC  data.  The  center  of  gravity  for  this  configuration  is  taken  as  3.45 
calibers.  The  static  margin  can  be  seen  to  be  about  0.6  calibers  at  Mach  number  2.00  and  to 
decrease  with  increasing  Mach  number  to  about  0.2  calibers  at  Mach  number  4.6. 

One  flight  datum,  from  a  limited  number  of  ballistic  range  stations,  seems  to  substantiate  these 
wind-tunnel  observations. 
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6.  THE  EFFECT  OF  WIND-TUNNEL  STING 

As  mentioned  in  Chapter  II,  the  NSWC  design  code  takes  into  account  the  effect  of  fin 
overhang  on  the  mutual  wing-body  influence  coefficient.  The  analysis  by  this  code  represents 
free  flight  conditions.  In  the  wind-tunnels,  the  models  are  held  in  place  by  sting  balances  that 
form  an  extension  to  the  body.  Since  the  diameter  of  the  balances  is  smaller  than  that  of  the 
body,  an  air  gap  is  formed  between  them  and  the  over  hanging  fins.  As  a  result,  the  analysis  of 
the  actual  situation  is  complicated  and  can  not  be  predicted  by  the  design  code.  Hence,  an  upper 
limit  for  this  effect  is  obtained  by  considering  a  configuration  with  the  body  extended  to  the 
trailing  edge  of  the  fins.  The  results  are  processed  and  only  the  effect  of  body  on  the  taii  is  taken 
into  account  (the  balance  does  not  include  the  loads  acting  on  the  sting).  The  upper  limit  esti¬ 
mates  for  the  increment  in  normal-force  curve  slope  and  the  backward  shift  in  center-of-pressure 
location,  due  to  the  balance,  are  shown  in  Figure  17.  The  shift  in  center-of-pressure  location 
reaches  0.2  calibers,  at  Mach  numbers  larger  than  4.0. 

If  the  results  of  the  wind-tunnel  tests  (Figure  16)  are  corrected,  to  exclude  the  effect  of  the 
sting  on  the  fin,  the  static  margin  reduces  to  0.5  calibers  at  Mach  number  of  2.0,  vanishes  at 
Mach  number  of  4.6  and  becomes  negative  (unstable)  at  higher  Mach  numbers.  These  values 
should  be  considered  conservative  because  the  correction  is  an  upper  limit  estimate  for  the 
difference  between  wind-tunnel  and  flight  configurations. 

7.  COMPARISON  WITH  ANALYSIS 

a.  The  ay  ial-force  coefficient  Axial-force  data  is  available  only  from  the  LaRC  tests.  The 
experimental  contribution  to  this  coefficient  is  obtained  by 

aca  -  ca  -  ca 

°BT  °SMOOTH 

The  data  for  the  extended  nose  configuration  was  used.  A  comparison  with  the  predictions, 
described  in  Chapter  IV,  Section  3,  is  presented  in  Figure  18  and  shows  a  good  agreement. 
Recall  that  the  comparison  is  valid  for  LaRC  wind-tunnel  conditions  and  that  the  full  scale 
axial-force  coefficient  due  to  the  buttress  threads  is  expected  to  be  much  larger. 

b.  Effect  of  the  buttress  threads  The  analysis  of  the  effect  of  the  buttress  threads  on  the 
normal-force  and  the  pitching-moment  coefficients  assumes  that  the  inner  sections  of  the  fins  are 
masked  by  the  boundary  layer  displacement  thickness,  associated  with  the  additive  axial-force 
caused  by  the  threads.  The  results  presented  in  Chapter  IV,  Section  3,  are  valid  for  free  flight 
conditions.  For  the  present  comparison,  the  analysis  was  repeated  for  the  LaRC  wind-tunnel 
conditions  with  the  pertinent  fin  efficiency  shewn  in  Figure  12.  The  slope  of  the  normal-force 
curve  is  practically  independent  of  the  roughness.  The  piedicted  forwaru  shift  in  the  center-of- 
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pressure  is  0.1  calibers  at  Mach  numbers  of  3.0.  reducing  to  0.05  calibers  at  a  Mach  number  of 
5.0.  These  changes  are  much  smaller,  especially  for  high  Mach  numbers,  than  those  predicted 
for  flight  conditions. 

The  experimentally  obtained  normal-force  curve  slope  is  shown  in  Figures  19a  and  20a  for 
the  open  and  closed  inlets,  respectively.  No  discemable  effect  of  the  buttress  threads  can  be 
detected  in  either  case.  The  effect  of  the  buttress  threads  on  the  center-of-pressure  location  is 
shown  in  Figures  19b  and  20b.  The  LnRC  results  show  a  0.1  caliber  improvement  in  static 
margin  due  to  the  buttress  threads,  in  the  case  of  the  open  inlet  and  no  effect  in  the  other  case. 
On  the  other  hand,  the  CRDEC  results  show  a  small  forward  shift  in  the  center-of-pressure. 
Taking  all  the  data,  together  with  the  uncertainties,  it  is  concluded  that  there  is  a  negligible  effect 
of  the  buttress  threads  on  the  normal-force  and  pitching  moment  characteristics  under  these 
wind-tunnel  conditions. 

c.  Ef.cct  of  ingested  inlet  flow  The  analytic  prediction  of  the  SFRJ  aerodynamic  characteris¬ 
tics  were  based  on  the  assumption  that  the  normal  force  for  the  open  inlet  configuration  is  exactly 
the  same  as  that  of  the  extended  nose  configuration.  This  in  turn  is  based  on  the  observation  that 
the  force  required  to  turn  the  ingested  air  is  the  same  as  the  slender  body  prediction  for  the 
extended  nose  tip. 

Figure  21  is  a  plot  of  the  data  discussed  in  the  previous  section,  showing  the  effect  of  the 
open  inlet  on  the  aerodynamic  characteristics.  It  is  observed  that  the  open  inlet,  with  air  flowing 
through  the  body,  has  a  higher  normal-force  curve  slope  than  that  of  the  configuration  with 
extended  nose  where  the  inlet  is  closed.  The  difference  is  about  0.49  at  Mach  number  2.0  and 
decreases  as  the  Mach  number  increases.  The  experimentally  obtained  difference  was  fitted  by 

ACm  "  *0.73  -  .124  M  ,  2.0  <  M  <  4.6  . 

The  center-of-pressure  location  with  open  inlet  is  0.34  calibers  more  forward  than  that  with 
the  extended  nose  tip.  This  shift  is  practically  independent  of  the  Mach  number. 


The  increment  in  pitching-moment  coefficient,  associated  with  the  increment  in  the  normal- 
torce  curve  slope  and  the  shift  in-center-of-pressure  is: 


12 


The  data  were  processed  and  the  following  curve  was  fitted  to  the  results: 

ACM  '*  -0.34  +  0.35  M  ,  2SMS4.6  . 

ma 

Within  the  range  of  interest,  both  ACM  and  ACm  are  positive,  indicating  that  the  center-of- 

ma  1,1  a 

pressure  of  the  contribution  due  to  the  open  inlet  is  upstream  of  the  lip  of  the  inlet.  This  implies 
that  the  additional  load  inc  udes  a  couple. 

d.  Discussion  The  predicted  aerodynamic  characteristics  of  the  60mm  SFRJ  projectile  having 
extended  fins  is  shown  in  Figure  21,  together  with  the  experimental  data  discussed  in  the  previ¬ 
ous  sub-section.  On  both  ends  of  the  Mach  number  range  the  experimentally  obtained  normnl- 
force  curve  slope  is  higher  than  that  predicted  by  1 1%  to  12%.  The  variation  of  the  experimental 
data  is  monotonic  and  does  not  show  a  jump,  as  predicted  by  the  NSWC  code. 

The  experimentally  obtained  center-of-pressure  location  of  the  configurations  having  the 
extended  nose  tip  are  0.25  to  0.40  calibers  more  aft  than  that  predicted.  Here,  too,  the  change  in 
location  between  low  and  high  super  onic  Mach  numbers  is  gradual.  The  gaps  in  normal-fore,: 
and  center-of-pressure  indicate  that  the  contribution  of  the  fins  to  static  stability  is  slightly 
underestimated  by  the  NSWC  code. 

The  center-of-pressure  location  of  t  ie  actual  configuration,  namely  the  one  with  the  open  inlet 
is  in  /ery  good  agreement  with  the  prediction  at  Mach  numbers  under  2.5  and  above  3.75.  In  the 
intermediate  range  the  experimental  curve  is  smooth.  It  appears  that  the  effects  associate  with 
the  open  inlet  offset  the  increased  contribution  of  the  tail. 

VI.  CONCLUSIONS 

An  hybnd  method  for  the  aerodynamic  analysis  of  the  SFRJ  projectiles  has  been  used  on  two 
si'in  stabil.zed  and  two  fin  stabilized  configurations.  The  method  uses  the  NSWC- 
Aeropredicrion  code  for  the  analysis  of  a  model  configuration,  having  an  extended  nose  and  a 
smooth  body.  Engineering  approximations  are  used  to  account  for  the  aerodynamic  effects  of  the 
inlet  and  the  buttress  threrdr 

The  NS  VC-Aeroprediction  code  was  selected  for  the  present  work  because  it  applies  empiri¬ 
cal  corrections  for  the  effect  of  fin  over-hang  on  tail- body  influence  coefficients. 

The  predicted  normal-force  curve  slope  of  the  40mm  SFRJ  is  18%-23%  lower  than  ex¬ 
perimentally  observed.  In  the  case  of  . the  75mm  SFRJ  the  difference  is  only  3.7%.  The  agree¬ 
ment  in  center  of-pressure  location  is  very  good  in  both  configurations. 
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The  hybrid  method  predicts  that  the  60mm  SFRJ,  with  the  original  Tins,  is  statically  unstable 
at  Mach  numbers  larger  than  2.0.  A  transonic  range  test  verified  this  prediction. 

The  analysis  predicts  that  the  60inm  SFRJ.  with  extended  fins,  has  a  static  margin  of  0.5 
calibers  at  low  supersonic  Mach  numbers.  That  diminishes  as  Mach  number  approaches  4.3,  and 
becomes  unstable  at  higher  Mach  numbers.  Wind-tunnel  tests  of  this  configuration  gave 
normal-force  curve  slope  which  is  1 1%  to  12%  larger  than  that  estimated.  The  experimentally 
obtained  center-of-pressure  location  is  0. 15  to  0.20  calibers  more  aft  than  the  analytical  value.  A 
study  of  the  results  show  that  the  contributions  of  the  open  cowl  to  the  normal-force  is  larger 
than  estimated  in  the  analysis  and  that  the  contribution  of  the  fin  set  is  slightly  larger  than  pre¬ 
dicted  by  the  NSWC-Aeroprediction  code. 

The  experimentally  obtained  contribution  of  the  buttress  threads  to  the  axial-force  coefficient 
is  in  good  agreement  with  the  predicted  value.  No  effect  of  the  threads  on  the  stability  of  the 
projectile  was  observed.  However,  the  expected  effect  in  the  wind-tunnel  flow  condition  is 
within  the  accuracy  of  the  measurements  and  is  much  smaller  than  that  predicted  for  free  flight 
conditions. 
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Figure  12.  Fin  efficiency. 
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Figure  13.  Aerodynamic  characteristics  of  the  60mm  SFRJ  projectile  having  orii  nai  tins. 


Figure  17.  The  effect  of  sting  balance  on  the  aerodynamic  characteristics. 


b.  Center-of-pressure. 


rv.  v  -u*  P.A  r\A 


60  mm  SFRJ 


8.0-i 


WftWWrWTXjrwnoWHfl  wiwr**v.nvrvrv*-i  ^iTS’rw>rw  tw  »■„  fytyrv  t.x-wtm  tkkvitu  t»  tukytl  tv  tvtvtw.wvtv 


NRSfl  LoRC  WIND  TUNNEL 
,  BUTTRESS  THREADS 


NASA  LoRC  WIND  TUNNEL 
*  BUTTRESS  THREADS 


i 

» 


'<  rr  ir»jr« jitw  %».  r*  r+  >ww  w»v ' 


K 

3 


5 

I 


Q 


Q/d°X 


ti 


o: 

Lk 

00 

1 

35 


o 

!/> 

u 

■a 

V. 

'C 

<u 

I 

-C 

o 

o 

'c 

3 

_§J 

£ 


8 


r-i 

K 

3 

W> 

LL 


K 

3 

VS 


-d 


Si 

i 


* 


44 


V^>*v 


% 
\ 
i 


V V V  V  V V  V  V  WV V  V ->-.*^>  ->«-, icym v-M-«k'  .  -  k'n^^in; 


REFERENCES 


!  1.  Vukelich,  S.R.,  Jenkins.  J.E..  "Missile  DATCOM:  Aerodynamic  Prediction  of  Conven¬ 

tional  Missiles  Using  Component  Build-Up  Techniques,"  AIAA  Paper  84-0388,  1984 

I 

k 

2  Mason,  L.A.,  Devan,  L.,  Moore,  F.G.,  McMillan,  D.,  "Aerodynamic  Design  Manual  for 
Tactical  Missiles,"  NSWC  TR  81-156,  1981. 

I 

i 

!  3.  Van  Dyke,  M.D.,  "First-  and  Second-Order  Theory  of  Supersonic  Flow  Past  Bodies  of 

!  Revolution,"  Journal  of  Aerospace  Science,  Vol.  18,  No.  3, 1951. 

|  4.  Syvertson,  C.A.,  Dennis,  D.H.,  "A  Second-Order  Shock-Expansion  Method  Applicable  to 

!  Bodies  of  Revolution  Near  Zero- Lift,"  NACA  Report  1328,  1957. 

|  5.  Nielsen,  J.  "Missile  Aerodynamics,"  McGraw-Hill  Book  Company  LTD.,  1960. 

1 

2  6.  Stoney,  A.,  "Collection  of  Zero-Lift  Drag  Data,"  NASA  TR  Report  100,  19nJ. 

t  7.  Hoerner,  S.,  "Fluid  Dynamic  Drag,"  Published  by  the  author,  1961. 

f  8.  Young,  A.D.,  Paterson,  J.H.,  "Aircraft  Excrescence  Drag,"  North  Atlantic  Treaty  Organiza- 

j  tion,  AGARD- AG-264,  198  i . 

; 

9.  Oskay,  V.,  'Transonic  Range  Test  of  the  40mm  SFRJ  Projectile,"  BRL  report  in  prepara¬ 
tion,  U.S.  Army  Ballistic  Research  Laboratory,  Maryland. 

10.  Olson,  D.N.  and  Wrede,  F.G.,  "Wind  Tunnel  Test  Results  on  a  34  Per  Cent  Scale  Model  of 
the  Fin  Stabilized  Solid  Fuel  Ramjet  Tubular  Projectile  at  Mach  Numbers  of  2  to  3.75," 

!  CRDEC  Technical  Report  to  be  published,  Aberdeen  Proving  Ground,  Maryland. 

11.  DATCOM-USAF  Stability  and  Control'Duta  Compendium,  Second  Edition,  1975. 

12.  Yalamanchili,  R.J.,  'Transonic  Range  Test  Results  of  the  75mm  SFRJ  Projectile,"  BRL 

!  report  in  preparation,  US  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 

j  Maryland,  1986. 

;  13.  ESDU  Data  Units  5.01.03.01  and  5.08.03.01,  Engineering  and  Scientific  Data  Units  Intema- 

|  tional.  Ltd.,  1983. 


45 


REFERENCES  (Continued) 


14.  Fenter,  F.W.,  "The  Turbulent  Boundary  Layer  on  Uniformly  Rough  Surfaces  at  Supersonic 
Speeds,"  Chance  Vought  Research  Center,  Report  No.  RF-E9R-2,  1959. 

15.  Voisinet,  R.L.P.,  "Influence  of  Roughness  and  Blowing  on  Compressible  Turbulent  Bound¬ 
ary  Layer  Flow,"  NSWC  Dahalgren,  Virginia,  TR-79-153,  1979. 

16.  Shapiro,  A.,  "The  Aerodynamics  and  Thermodynamics  of  Compressible  Fluid  Flow,"  Vol. 
II,  Chapter  27,  The  Ronald  Press  Company,  1954 

17.  Yalamanchili,  R.J.,  Wood,  R.,  Test  results  from  NASA  Unitary  wind-tunnel  on  SFRJ  to  be 
published,  Aberdeen  Proving  Ground,  Maryland,  October  1987. 

18.  Olson,  D.N.,  Wrede,  F.,  Test  results  from  USA  CRDEC  supersonic  wind-tunnel  to  be 
published,  Aberdeen  Proving  Ground,  Maryland,  October  1987. 


vuim  b*njwwuurnjf»rii.ir*  trm 


wn  mt krwn wrvrrvvAJVvv\*ru^n  r/lw  mu'*  *  *:*■¥  m* wur wm w\  wntimw w.v\v  -*a ttv*  \nrvr  !j*  it*  v *  ■  r  v km s 


LIST  OF  SYMBOLS 


AR  aspect  ratio 

B  VM2  -1 

BT  buttress  threads 

c  chord 


root  chord 


CA  axial-force  coefficient,  based  on  reference  area 
Cpg  drag  coefficient  of  an  excrescence,  based  on  frontal  area 


CN  normal-force  coefficient 


CN(x  normal-force  curve  slope 


CM  pitching-moment  coefficient 


diameter  of  center  body 


Dj  diameter  of  inlet 


h  height  of  a  protuberance 

H  boundary  layer  shape  factor 

KT(B)  influence  coefficient-tail  in  presence  of  body 

Kg(T)  influence  coefficient-body  due  to  tail 


/,  length  of  extension  of  the  nose 


M 


Mach  number 
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LIST  OF  SYMBOLS  (Continued) 


S  area  of  exposed  pair  of  fins 

S  area  masked  by  the  boundary  layer 

SR  reference  area  jiD2/4 

i 

Xcp  center-of-pressure  measured  from  inlet  lip 

Greek  Symbols 

angle  of  attack 
boundary-layer  thickness 
boundary-layer  displacement  thickness 
boundary-layer  momentum  thickness 

tail  efficiency  factor 

Subscripts 
b  body 

e  extended  nose 

f  fin  assembly 

i  inlet 

o  open  inlet 

t  effective  flare 
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LIST  OF  SYMBOLS  (Continued) 


>uperscripts 


series  of  protuberances 


conected 


difference  between  open  inlet  and  extended  nose  tip 
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